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ABSTRACT 



> 

t:;;^- , We discuss the evolution of cosmological gamma-ray burst remnants, consisting of 

I the cooling and expanding fireball ejecta together with any swept-up external matter, 

' after the gamma-ray event. We show that significant optical emission is predicted 

■ which should be measurable for timescales of hours after the event, and in some cases 

G^ ■ 

radio emission may be expected days to weeks after the event. The flux at optical, 
Q^. X-ray and other long wavelengths decays as a power of time, and the initial value of 

the flux or magnitude, as well as the value of the time-decay exponent, should help to 
distinguish between possible types of dissipative fireball models. 



Subject headings: gamma rays: bursts 



1. Introduction 

Gamma-Ray Bursts (GRB) must leave behind remnants, hereafter referred to as GRBR. 
In general, GRB at any plausible distance should result in relativistically expanding fireballs 



(e.g. Meszaros, 1995 for a review), and the 7-ray emission most probably arises after the fireball 
becomes optically thin, in shocks occuring either beacuse the ejecta run into an external medium, 
or because internal shocks occur in a relativistic wind. In the first type (a) of GRB model, the 
initial energy input is impulsive and the relativistic ejecta begin to decelerate when they have 
swept up external mass that is a fraction ~ of the ejecta mass, leading to an external and 



a reverse shock that radiate a large fraction of the initial total energy (the GRB event, Rees fc 
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Meszaros, 1992, Meszaros & Rees, 1993a, Katz, 19*9^ , pari, Narayan fc Piran, 1996| ). However, 



after this initial burst the ejecta should still continue to expand, sweeping up an increasing amount 
of matter and slowing down. As pointed out by Paczyiiski &: Rhoads, 199^ , the evolution of such 



an "old" cooling fireball should resemble the Van der Laan model of expanding radio sources, and 
they estimated the possibility of detecting the late radio emission from such objects. In the second 
type (b) of GRB model ( Rees fc Meszaros, 1994 , [Paczyhski &: Xu, 1994|) , the initial energy input 



continues over a period of time and the resulting relativistic fireball wind produces a GRB due 
to internal shocks in the wind itself. These objects too should leave behind a cooling remnant, 
although the dynamics of the late evolution are expected to be different from those of type (a) 
because the burst occurs at significantly smaller radii and the effects of the external matter are 
negligible for some time afterwards so the dynamics are different. 

In this paper we investigate the dynamical evolution of the GRBR following the GRB 
event, both for the impulsive external shock and the wind internal shock models. We assume 
that the GRB are at cosmological distances (similar models with somewhat different radiation 
characteristics can be calculated for galactic halo models too). Both types of cosmological GRBR 
produce, in addition to decaying 7-rays, significant amounts of softer radiation, mostly X-rays 
and optical, but in some cases also radio. Depending on the distance to the source, the duration 
of the main GRB event and the physics of the burst model, this radiation can be in many cases 
detectable with appropriate instrumentation. Based on practical considerations of sensitivity, 
telescope availability and fast response time, the best chances for detection at energies other than 
gamma-rays may be at optical wavelengths. 



2. Dynamical Evolution of GRBR 

A source at distance D expanding isotropically with a relativistic bulk Lorentz factor F, which 
in its own frame has a comoving specific intensity I^^ at comoving frequencu z/^, will produce as 
a function of observer time t an observer-frame flux at frequency Vm 

where I'm = ^^'m where the observed transverse apparent size ctV is used ( [Rees, 1966| ). In 
the initial phase of expansion of the flreball the internal energy density drops adiabatically and 
most of the energy is locked up in the expansion kinetic energy, so that I^^ is negligible until 
after shocks dissipate some or all of the bulk kinetic energy in the GRB event. At this point, 
the GRB spectrum is generally considered to be due to synchrotron or inverse Compton (IC) 
radiation from power-law electrons or pairs accelerated in the shocks. A break in the spectrum 
appears at photon energies associated either with synchrotron photons whose index changes below 
and above the minimum synchrotron frequency corresponding to electrons with 7mm) or else with 



IC photons corresponding to the same break but shifted up in energy by (Meszaros, Rees &: 



Papathanassiou, 1994, henceforth MRP94). If electrons have an energy distribution oc 7 above 
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7mm) the uFi, photon spectrum below the break is oc (or i/^/^, but we shall take oc z^"^ as a 
generic observed spectrum), and is flat above that. The GRB spectral breaks which are detected 



by BATSE generally in the 0.1-2 MeV region, e.g. Band, etal., 1993 (although they could have 



a much broader distribution, Piran &: Narayan, 1996| ), could be due to the synchrotron break 



(e.g. the PC model in RMP94, or [Katz, 1994| , |Sari, Narayan k Piran, 1996| ). Alternatively, a 



synchrotron spectral break might appear at significantly lower frequencies, and the similarly 
shaped IC scattered spectrum could provide the MeV break (MRP94, models FC and TC). A 
key quantity therefore, both for the GRB and the ensuing GRBR, is the comoving synchrotron 
intensity at the comoving synchroton peak frcQUcncy ~ W^B'-fl^i^Rz, where B' (G) is the 
comoving magnetic field strength, 

Il oc ^'^^ 7"""'^^ oc n'.B'AR' . (2) 

imin 

The GRB from impulsive fireball external shocks are expected to produce simultaneous X-ray 
and optical flashes ( Meszaros &i Rees, 1993b| , MRP94) whose duration is similar to that of the 



gamma-ray burst. Continuous wind flreball shocks also produce flashes of softer radiation which 
are contemporaneous with the GRB flash ( Papathanassiou fc Meszaros, 1996| ). The relative 



amount of contemporaneous soft emission depends on the type of shocks, the particle acceleration 
process and the strength of the magnetic field in the acceleration region. For the GRBR, the shock 
physics is also a major input in determining the comoving intensity I^^ (equation [l|), and thus 
the relative amounts of soft emission expected after the main gamma-ray and contemporaneous 
X-ray/optical flash is over. To a very large extent, however, the temporal behavior of the GRBR 
will be dominated by the evolution of F and some of the related dynamic quantities, which are 
different depending on the initial energy input regime giving rise to the fireball. 

a) In the impulsive model, the flreball initially expands at its saturation bulk Lorentz factor 
F ~ 7/ = (Eo/MoC^) ~ constant, and upon interaction with an external medium it develops 
a blast wave moving into the external medium and a reverse wave moving into the ejecta, 
which radiate away a substantial fraction of the initial energy Eq to give the GRB. This occurs 
at a relatively large radius r^ec ~ 10^^(-E'5i/next)^''^??3 ^''^cm over an observer-frame timescale 
^7 ~ tdec ~ i^decC^^ '^ ~ 1 (-£'51 /"-exi ) ^''^??3 ^^^s, which is in the range 1 — lO^s for r/3 ~ 1 — 10~^ 
(this is the observer-frame burst duration, for an impulsive type of initial energy input occuring 
over a timescale shorter than tdec)- During the GRB the radiative timescales are short enough 
to ensure high radiative efficiency, but the initial rapid cooling is generally sufficient to ensure 
that the subsequent evolution of the remnant will be adiabatic. For such impulsive fireballs, if 
the external medium is approximately homogeneous, the GRBR will continue to evolve afterwards 
with 

F oc oc ; r ^ ctT^ o, t'/\ (3) 

Here F is the Lorentz factor of the contact discontinuity (CD) between ejecta and external 
medium, and essentially also that of the external shock, r is the distance advanced by the CD 
along the line of sight (longitudinal size), and t is observer- frame time. The comoving expansion 
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(adiabatic cooling) time is tg^ ~ rc ^ oc t^/^, and the comovinff radial extent of the shell of 
fireball and swept-up matter is AR' ~ ct'^^ oc t^/^. 

b) In the continuous input (wind) model, before the GRB the wind is expanding at its saturation 
average bulk Lorentz factor F ~ = {L/MqC^) ~ constant, and the primary energy input remains 
continuous (rather than being impulsive) over a wind duration timescale tw which characterizes 
the burst duration in the observer frame. Variations on a timescale t^ar < tw of the bulk Lorentz 
factor Ar ^ T, e.g. due to corresponding L or M variations, lead to internal shocks at relatively 
smaller radii r^j^ ~ ctvar^'^ ~ 3 x 10^^t„ar,-3??2 cm, which randomize the relative kinetic energy of 
different shells of the wind and radiate away a fraction of order unity of the total wind kinetic 
energy (the GRB event; e.g. |Rees &: Meszaros, 1994 , Paczyiiski Sz Xu, 1994 , Waxman fc 



Piran, 1994). During the shocks (for t ^ t^) F does not change very much and may be taken 
as approximately constant. After the shocks cease {t > t^), the wind remains at F ~ ary ~ 
constant (with q ^ 1 ~ constant) as long as deceleration by an external medium does not come 
into play. (If and when this occurs, an external shock should develop, leading to a slow motion 
version of an impulsive GRB; the remnant thereafter behaves similarly to to case (a) above). In 
the range rdis ^ r ^ r^ec (where r^ec ~ 2 x lO^''^(ii^5i/ne^t)^/^0~^/^r/2 '^^^ cm is reached after a time 
idee ~ {E^i /next)^^^G~'^^^r]2 ^^^s) the wind and the remnant evolve with 

F ~ a?7 constant ; r ~ ctF^ act , (4) 

where as before r is the distance advanced by the CD (the front of the wind) in observer time t. 
During the wind regime and for t < t^, while internal shocks continue producing the GRB light 
curve, the comoving density evolves as n' oc cx and any comoving wind (transverse) field 
as B' oc oc t~^. However, for r ^ rj^p ~ ct^F^ ~ 3 x 10^'^twr]2 cm, or for times ^ t ^ tdecj 
the wind has stopped blowing while the flow can still have (temporarily at least) F ~ constant 
and r (X t. The flow is now in an impulsive regime (observer-frame durations are ^ tw), and the 
gas expands isotropically in its own rest frame. The very rapid cooling during the GRB again 
ensures that the subsequent evolution of the smoothed-out remnant wind is adiabatic. In this 
GRBR regime then n'^ oc oc B' oc r~'^ oc and particle random Lorentz factors cool 
adiabatically as 7 oc n'^^^ oc oc t^^. 



3. GRBR Spectra from Impulsive Fireballs 

One can think of three main sub-cases which share the impulsive external shock dynamics 
described in §2 (a), but which, depending on the type of shock physics involved lead to different 
GRBR spectral evolution regimes. 

(al) For the simplest impulsive model, only the forward blast wave radiates efficiently (e.g. the 
reverse shock might be an inefficient particle accelerator or an inefficient radiator; this might be 
the case if the reverse shock were very weak, perhaps due to a very high Alfven sound speed). 
This corresponds to the PC (piston) impulsive GRB model described in RMP94, e.g. their 
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Figure Ic. As a specific example, this model has a 7-ray 

erg cm~^ near the BATSE threshold, where the total energy is 10^^ £"51 erg channeled into an 
angle lO~^0_i radians at a luminosity distance 10^^D28 cm (roughly 3 Gpc; the same object 
at 300 Mpc would give one of the brighter bursts). This GRB spectrum satisfies amply the 
X-ray paucity constraint {Lx <C lO^^L-y), being due to synchrotron radiation that peaks at 
Ebr ~ 0.5r/35'7^.„ MeV, with a i^F^ slope near 1 below Ei,r and near above it. The behavior 
below the break extends during the GRB at least down to optical frequencies, but becomes 
self-absorbed well above the radio range. This spectrum would be expected from electrons 
with an energy power law index p ~ 3 above 7mm ~ kT ~ 10^7/3^2 ■ Here rj = 10^% is the 
final coasting bulk Lorentz factor, i?' ~ 10 G is the turbulently generated comoving magnetic 
field strenght and 7mm is the minimum post-shock electron random Lorentz factor, assuming 
that electrons and protons achieve a fraction Kme/nip of their equipartition energy (for details 
see MRP94). In the GRB event, the fireball loses on the order of its initial kinetic (total) 
energy, and we can assume adiabatic conditions for the GRBR evolution, i.e. the comoving 
cooling time rapidly becomes longer than the comoving expansion time. The comoving intensity 
dll) from the external matter that produced the initial GRB and which subsequently evolves 
adiabatically is I^^ oc t^^/^, since B' oc V'~'^^^, n'^ oc V'^^ oc t^^/®, 7 oc r^'^ oc t^-*^/^. Therefore 
Fy^ OC t^T^Il^ oc t~'^/^ starts to decrease in time immediately after the GRB for this initially 
shocked shell of external matter, its peak frequency dropping as Vm oc TB'^^ oc t^^^/*^, giving 
Fy^ oc z/m^^. This adiabatically cooling leftover component from the GRB is, however, weaker 
than that due to newly shocked external matter, as the ejecta continues to advance at a steadily 
decreasing velocity into the external matter. For this newly shocked matter B' , rig, 7 are all 
three independently oc F oc so /^^ oc As a result, from (|l|) the observed flux at the 

break evolves initially as Fy^^^ oc t*^ ~ constant, so the flux below the break remains constant, 
Fy (X u ~ constant as long as u < Vm, but the latter decreases in time as i^m oc we 
take approximately the gamma-ray band to be at 10^^ Hz and the optical band at 10"^^ Hz, i^m 
drops below the optical band at topt ~ 3 x lO^t^, where t-y ~ tdec is the duration of the gamma-ray 
flash defined below equation (||). Thus, since the flux is oc above fm the detected optical flux 
from an impulsive (al) type GRBR is Fopt ~ 10~^^ K = constant for t ^ topt ~ 3 x lO'^t^, 
and Fopt ~ W^^'^K D2^{t/topt)~^^'^ for t ^ topt, where flux units are erg cm^^ s^^ Hz^-*^, D is 
the luminosity distance and K = E^iOzf is the burst energy normalization . The corresponding 
V-band magnitudes are m„ = — 2.5 log(Fopt/4 x lO"^'^) or 

~ 19 -2.5 log A" + 5 log 1)28 , for t ^ 3 x lO^t^, , 

(5) 

~ 19-2.51ogK + 51ogi:»28 + 3.751og(t/topt) , for t ^ 3 x 10^^ , 

The optical brightness is constant at the level given by the first line of equation (^) for about an 
hour (since for the fiducial parameters here ~ 1 s). Note that at 300 Mpc (-D28 = 10~^) the first 
line of (P) would be ~ 14 + • • •. In this model, the self-absorption frequency (initially at ~ 10^^ 
Hz) evolves at first oc slower than t'm, the two becoming equal at tm,afe 4 x lO^t-^ where 
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t'ab ~ 3 X 10^^ Hz. Thereafter, the frequency Uab oc t ^^/^^ and the peak of the spectrum separating 

the optically thick and thin regimes evolves as Fi,^^^ oc i-']^^^'^ as long as the GRBR still expands 
relativistically (it reaches F 1 after t^r ^ lO^tj, about three years for '-^ 1 s), and F^^^ oc v^j^ 
afterwards. As a result, the radio flux in the 100 MHz band is neglig ible, Fr ^ IQ-'^KD^inJy. 
(This refers to the usual incoherent self-absorbed synchrotron flux; the black-body upper limits 
used for this estimate could of course be violated by coherent emission mechanisms). 

a2) A slightly more involved impulsive model considers both the reverse and the forward shock to 
be efficient radiators. An example of this occurs if for instance the ejecta has frozen-in magnetic 
fields, sufficiently strong to ensure radiation but not enough to dominate the dynamics (the FC 
model of RMP94, their Fig. la). At t-y ~ tdec the reverse shock has achieved transrclativistic 
speed (fV ~ 2 respect to the ejecta, with k = 10^, B' rsj 0.3G and F ~ r/ ~ 10^ in this 
numerical example). The ejecta electrons have 7mm 10^, giving an observed synchrotron peak 
at ~ 3 X 10^"^ ~ 10^^ Hz with fluence 5"^ - lO'^KD^i erg cm ex. above Vm (or flux 

oc i^-^), where K = E^iOZh then between 10^^ and 3 x 10^° it has Su ~ (flux ~ i/°), 
and Si/ '~ erg cm-^ oc (F^ oc ^) above 3 x 10^° Hz, the latter two corresponding 

to IC-scattcred synchrotron by reverse shock electrons. We consider the evolution of the reverse 
shock (ejecta) region only, since it has stronger XR and O contributions than the forward blast 
wave. (The forward blast wave remnant might be naively expected to behave as the al case above, 
but being an IC component its contribution drops in time by two powers of 7 faster than the 
synchrotron component). The shocked reverse gas is in pressure equilibrium with the gas forward 
of the CD, so it also moves with F ~ t~^^^ and r oc t^^^ as does the CD. The ejecta frozen-in field 
behaves as B' oc oc t~^^'^. The comoving pressure on either side of the side of the CD evolves 
oc so the reverse gas, even though traversed by reflected pressure waves, cools in its rest 
frame with T' oc r~^/^ and the comoving density in the ejecta evolves as n' oc r~^^^ oc t-9/20, the 
total number of particles in the ejecta remaining constant. The ejecta comoving radial width is 
AR' oc r~^/^ and the ejecta column density is n'.AR' oc oc t~^/'^. As the ejecta cool, any fresh 
acceleration of ejecta electrons can at most give 7 oc nTp oc f^^^^, so the GRBR spectrum will 
be dominated by the adiabatic cooling of the electron energy acquired in the initial GRB shock 
heating, 7 oc r^^ oc t^^^'^. The reverse peak i^m ^ 10"'^^(t/t^)^"'^^/^ Hz is produced by a comoving 
peak intensity oc n'^B'AR' ~ lO'^it/t^)-'^. Therefore F^^ oc t^F^/^^ ~ W-^'^^KD^iit/t^y^/^ 
evolves proportional to vU^^- The self-absorption frequency initially is Vab ^ 10^^(1/^)^^/^ 
Hz, and the absorption and synchrotron peak frequencies become equal at v^m^ab ^ 10^^ Hz 
after a time t^^ab ~ lO^t^. Thereafter the peak of the self-absorbed flux occurs at a frequency 
i^ab ^ 10"(t/t„t,a6)"^'^/^^ where F^^^ - 3 x W'^^iu/uabf^^^'^ ■ The optical flux near 10^^ Hz is 
Fopt ^ 10^^'^(t/t^,)^^/^ for t > t^, both before and after tm^ab- (In principle one might have 
expected somewhat later a phase of Fgpt ~ constant due to the oc IC component coming into 
the optical range, but as mentioned above, the IC energy losses drop off in time much faster than 
the synchrotron losses, so this should not be important). The visual magnitude of the GRBR from 
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this reverse shock with frozen-in magnetic fields is therefore 

m^, ~9-2.51ogK + 51ogL>28 + (45/8)log(t/^) , for t ^ . (6) 

The incoherent flux in the 100 MHz radio band is maximum when the self-absorption peak passes 
through that frequency at the time tn ~ lO^t^, giving Fy^^j^ax ~ IKD^in^y, which at 300 Mpc 
would be ~ O.IKD2Y niJy, still very low. Before the maximum of the radio flux (between t^^ab 
and t{i) the radio flux grows as F^^^ oc t^^/-*^^, and after it decays as F,y^ ~ 

a3) A somewhat similar impulsive model also considers both reverse and forward shock emission, 
but differs from (a2) in the origin of the magnetic field, which is here assumed on both sides of 
the contact discontinuity to be built up by turbulent motions to some fraction of the thermal 
proton energy density. The GRB spectrum corresponds to the TC model of RMP94 (their Fig. 
lb). In this example the field energy is a fraction A ~ 10~^ of equipartition, and all initial 
physical parameters (as well as the initial spectrum) are the same as in case a2) above; the 
optical peak is synchrotron from the reverse shock. However, the magnetic field evolves in time 
differently, S'^/Stt maintaining itself roughly at a constant fraction A of the thermal proton 
energy density kn'T' oc oc t~^/^ , so B' (X in the reverse shocked ejecta. Retracing 

the steps in a2) we now have I^^ oc t^"^^^, Vm t~^f^, Fy^ oc t~^/^, or Fy^ oc initially. 
Also Uab oc and becomes coincident with at Vm,ab ~ 10^°'^ Hz at the time 

tm,ab ~ 3 X lO^t^, where the peak value is Fy^^^ ~ 3 x 10^^^. After that time the peak of the 
self-absorbed spectrum moves down as Fy^^^ oc v^^^^^ , with oc ^-53/56^ ^pj^g optical flux around 
lO^s Hz is Fopt ~ 10-23 /i:£)-^(t/t^)-2 

erg cm ^ s ^ Hz ^ both before and after tmafej so the visual 

magnitude of the GRBR is 

~9-2.51og/C + 51ogL>28 + 51og(t/^) , for t ^ . (7) 

The self-absorption peak comes into the 100 MHz radio range at t/j ~ 3 x lO^tm^af, ~ lO^t^ 
(a few weeks after the GRB, for fiducial parameters with ~ Is), and the peak radio flux is 
Fy^^max ^ lO-^KD^g^ niJy, or more interestingly, Fy^^^nax ~ l-f^^27^ mJy for a GRBR at 300 
Mpc. Between tm,ab and the radio flux grows as Fy^ oc t^^/-*^^, and decays as Fy^ ~ t"^ after Ir. 



4. GRBR Spectra from Wind Fireball Models 



For wind models, GRB spectra have been discussed, e.g. by Meszaros fc Rees, 1994 and 



Thompson, 1994 , [Thompson, 1996 . The observed fluxes from the resulting GRBR will be given as 



before by equations (1) and (2), but there are signiflcant differences from the impulsive case due to 
the smaller radii at which shocks occur, the different evolution dynamics (§2, b) and the different 
post-shock behavior of the physical quantities. Below we consider two particular examples of wind 
internal shock models. 

bl) The simplest wind GRB spectrum would be one where the break around lO^'^Hz ~ 0.5MeV is 
due to the synchrotron peak, e.g., a case resembling the (al) impulsive spectra. We take as a specific 
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example ( Papathanassiou fc Meszaros, 199q ), with T] = 10^, = lO^S, t^ar ~ 2 X IQ-^s, 



and at r^is ~ 6 x 10^"^ tyar -2112^''^ the turbulently generated field is B' ~ 10^ G and k ~ 10^ 
(internal shocks have Tj-ei ~ 1, so 7mm ~ 10^ and Um ~ ~ lO^'^ Hz. The total energy into 
9 ~ lO~^0_i is again E = lO^^-Esi erg spread over tu, = 10^t^^2S. The GRB spectrum has the same 
shape and fluence as the impulsive case (al), but the flux level is 10~^ because the duration is 10^ 
longer in this example. We consider the GRBR regime for r ^ rj^p ~ 3 x 10^^t^2f/2^iii. However, 
in this type of model the comoving t'^y^^^ ~ 5 x 10^''i3g"^73^^jj„ <^ tg^,, and the remnant undergoes 
essentially instantaneous cooling as soon as shocks stop heating the gas. This is characteristic of a 
wind model with a GRB break at 0.5 MeV ascribed to a synchrotron peak, since this requires high 
B' and 7. So in these cases we expect no significant emission from the GRBR after the gamma-ray 
flash, although there will be (as in al) some contemporaneous optical and X-ray emission during 
the GRB itself, for t <tw 

b2) Consider now GRB from wind internal shocks with a synchrotron peak at optical frequencies 
and an IC peak around 0.5 MeV. As a specific example ( Papathanassiou fc: Meszaros, 1996| ) we 



take a model with t^, ~ lO^s, t^ar ~ 2 x 10~^s, k = 10^, r/ = 10^, and at r^is ~ 6 x lO^^tt, 
we assume that magnetic fields build up by turbulent motions to a fraction A ~ 10~^'^ of the 
proton thermal energy, or ~ 3 G, which leads then to a GRB synchrotron optical peak at 
rjv'^ ^ 10^^ Hz, and an IC peak near lO^'^Hz ~ 0.5 MeV. The GRB fluence spectrum vF^t^u is the 
same as for the impulsive models b2and b3, and the flux has the same shape but is a factor 10~^ 
lower since here ~ 10^ s. In this case, the comoving synchrotron timescale t'^y ^ 5 x 10'^ ^ t'^^, 
so the adiabatic approximation will be valid in the remnant, which retains a significant thermal 
energy content that is radiated away during the later expansion. During the GRB the initial 
optical peak at ~ lO^^Hz has flux F,^^ r-j 10~^^t~2-^-^28^ Hz~^. The comoving 

intensity is, from equation (1), /^^ = F^^D"^ / (c^t^T^) ~ 10~'^ K D2gt~2V2^ ^'^d the self-absorption 
frequency is T{I,y'C^ ~ lO^^Hz. For the time scaling we use the properties of the 
wind in the impulsive F ~ constant phase (Q), so the comoving synchrotron intensity scales as 
II^ ~ n'^AR'B' (X t~^, the synchrotron peak frequency as Um ~ Fi?'7^ cx t~'^, and the observed flux 
at the synchrotron peak as F^^ ~ (P't^D'^T^Il^ cx so that the synchrotron peak moves down 
with peak frequency as F^^ cx vll"^ . The self-absorption frequency is Vah oc F(/^/7)^/^ cx 
and the synchrotron maximum and the self-absorption frequencies coincide at tm,ab ~ 10^ ^^tw at 
^m,ab 2 X 10^" Hz, where the flux is Fy^^^^ ~ 10"^''. After that the self-absorbed peak of the 
spectrum moves down in frequency as Fnu^b ^ ^tib' ■ '^^^ optical flux as a function of time is given 
by Fopt ~ Fy^{vm/IQ^^) ~ IQ^'^^KD^iit/t^)-^. The corresponding GRBR visual magnitude is, 
as a function of time, 

m„ ~ 14-2.51ogi^ + 51ogi?28 + 151og(Vt«,) , fort > , (8) 

where K = E^i9zit^\- The initial magnitude (14, in this example) is constant for t ^ = 10^ s 
during the GRB, and it drops extremely fast in the GRBR phase afterwards. AT 300 Mpc, the 
initial magnitude is my ~ 9, dropping to 24 after t ~ lOt^ ~ 10^ s for these parameters. The radio 
flux in this model is very low, Fr ~ ^.'^>KD2i^J.iy, or 0.03i^D^8^ mJy at 300 Mpc, growing as t 
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before tji ^ 3 x 10^^,^, ~ 3 x lO^s in this example, and dropping as t ^ afterwards. 



5. Discussion and Observational Prospects 

Gamma-ray burst remnants, or GRBR, should according to our calculations leave behind 
an afterglow at wavelengths longer than 7-rays, in particular at X-ray, optical, and in some 
cases also in the radio bands. The calculations discussed above apply both to spherically 
expanding configurations and to jet configurations, as long as the jet opening angle 6 ^ F^^. 
The adiabatic approximation made for the evolution of the remnant should be generally valid for 
times significantly longer than the gamma-ray burst duration. A wide variety of models, involving 
several unknown parameters, are compatible with the gamma-ray data. Observations in other 
bands offer the chance to narrow the range of options and refine existig models. The numerical 
examples discussed are illustrative of the possibilities, and should give a reasonable idea of the 
range of values that might be expected. There is greater reason for confidence in the scaling laws 
discussed, and in particular in the time dependence laws of the fiuxes in different regimes for the 
various models, as these are based on simple physical arguments. 

Omnidirectional, or at any rate large field-of-view X-ray detectors in space (such as HETE, 
Ricker, 1992| ) may be able to detect the simultaneous X-ray emission predicted for GRB (e.g. 



MRP94 jPapathanassiou fc Meszaros, 1996 ) as well as, for the brighter bursts, possibly the X-ray 



afterglow implied by our models discussed here, although surface areas higher than HETE's may 
be required to follow the remnant evolution. 

The synchrotron radio fluxes from GRBR are expected to be very small, reaching their 
maximum value on timescales of weeks after the GRB outburst, which in the optimal cases (e.g. 
§2, b3) may be at most in the mJy range. The radio fiux is expected to change before and after 
the maximum as power laws of the time which are characteristic of the models. Nonetheless, radio 
searches, even for very short timescale flashes, are worthwhile because we cannot rule out coherent 
emission behind relativistic shocks, which would of course permit brightness temperatures far 
higher than the usual self-absorption limit, (indeed, the intraday radio variations in AGN [[Wagner 



|1996[ may exemplify coherent emission from shocks whose properties resemble those expected in 
GRBRs). 

The optical detection of the GRB event itself is within the range of capabilities of modest size 
telescopes, but the problem is one of field of view, as typical BACODINE error boxes supplied 
within a minute of the GRB event are at least several degrees, and sometimes as much as ten 
degrees wide, and any significant improvement in location takes at least days. Thus detection of 
optical emission from the burst itself (rarely longer than minutes) is difficult. However the GRBR 
optical emission decays on longer timescales, typically hours, and for the initial magnitude levels 
implied by several of the models even at Gpc distances (e.g. 9 — 14, see §3 a2, a3, §4 b2), 

meter class telescopes equipped with CCD detectors may be able to cover a several square degree 
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wide field by doing rasters of ten minute observations. The advantage of such observations, if they 
can be repeated more than once over the same field, is that they would allow a determination of 
the optical flux time-decay exponent, which can help discriminate among models. 

We are grateful to NASA NAG 5-2857 and the Royal Society for support. 
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